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Abstract 

Future neutrino oscillation experiments will lead to precision measurements of neutrino mass splittings 
and mixings. The flavour structure of the lepton sector will therefore at some point become better known 
than that of the quark sector. This article discusses the potential of future oscillation experiments on 
the basis of detailed simulations with an emphasis on experiments which can be done in about ten years. 
In addition, some theoretical implications for neutrino mass models will be briefly discussed. 
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1 Introduction 

The observation of atmospheric neutrino oscillations by the SuperKamiokande experiment [J triggered a 
remarkable discovery phase. The initial evidence turned into a solid proof of neutrino flavour conversions 
as well as of the L/E dependence as required by oscillations. The solar neutrino problem has also been 
resolved in the last years. The Gallex experiment 0 detected initially a rate effect which implied flavour 
conversion on the basis of solar models. The SNO experiment proved then model independent neutrino 
flavour transitions [3 El The initially allowed parameter islands were cleaned up by the KamLAND ex¬ 
periment, which demonstrated finally with reactor anti-neutrinos [3 El that the so-called LMA-solution is 
correct. Altogether the existing experimental results fit now very nicely into a picture with three massive 
neutrinos, which corresponds to the simplest scenario for three generations. The only exemption is the dis¬ 
puted LSND result 0, which would have far reaching consequences if it were confirmed, but this possibility 
will be ignored here. The oscillations of three neutrino generations involve then two mass-squared differences 
A mf 2 — an d Am| 3 ~ A, three mixing angles, #12, $23, and $13, and a CP-violating phase 5 . 

Atmospheric neutrino data [H] and the first results from the K2K long-baseline accelerator experiment 3 
determine Ato| 3 = (2.2+q|) x 10 -3 eV 2 and 6*23 ~ 45° (HUS], whereas solar neutrino data mm , com¬ 
bined with the results from the KamLAND reactor experiment 0 lead to Am 2 2 = (8.2+^3) x 10 5 eV 2 
and tan 2 $12 = 0.391 qq4 0. The results can now approximatively be summarized by two independent two 
flavour oscillations where the solar and atmospheric Am 2 values are roughly now known. 

The key parameter for genuine three flavour effects is the mixing angle $i 3 which is so far only known to be 
small from the CHOOZ [T7H i.‘ji and Palo Verde |2| experiments. The current bound for 9 \3 depends on the 
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value of the atmospheric mass squared difference and it gets rather weak for Am^ < 2 x 10 ~ 3 eV 2 . However, 
in that region an additional constraint on 6*13 from global solar neutrino data becomes important IlSj . At 
the current best fit value of Am^ = 2.2 x 10 -3 eV 2 we have at 3cr the bound sin 2 $13 < 0.041 UJ. There is 
no reason why $13 should vanish and one should expect therefore $13 to be finite. 

One might think that neutrino oscillations are in future less interesting, since it will lead only to parameter 
improvements of the leading 2x2 oscillations and maybe a finite parameter value of $ 13 . However, such a 
view misses completely the fact that the neutrino sector is, unlike the quark sector, not obscured by hadronic 
uncertainties. The precision to which the underlying flavour information is determined will therefore only 
be limited by the ultimate experimental precision. If high precision measurements are possible, then they 
will be very sensitive tests of flavour models and related topics, like the unitarity of three flavours. Genuine 
three flavour oscillation effects occur only for a finite value of $13 and establishing a finite value of $13 is 
therefore one of the next milestones in neutrino physics. Leptonic CP violation is another three flavour 
effect which can only be tested if $13 is finite. The usual see-saw scenario includes besides 6 in addition 
two further Majorana CP phases in the light neutrino sector, as well as other CP phases in the heavy 
Majorana sector, which are involved in leptogenesis. In general the heavy and light CP phases are not 
connected, but most flavour models create relations between these two sectors, relating thus low energy 
leptonic CP violation to leptogenesis and mass models. Precision measurements of neutrino oscillations test 
therefore very interesting questions of particle physics which are connected to the origin of flavour and to 
phenomenological consequences of flavour. There is thus a very strong motivation to establish first in the 
next generation of experiments a finite value of 0\ 3 in order to aim in the long run at a measurement of 

leptonic CP violation ca El ESI ESI EDI- 

2 Three neutrino oscillation in matter 

An effective two flavour treatment is insufficient for future oscillation experiments and matter effects must 
be included in addition. The generalization of the oscillation formulae in vacuum to N neutrinos leads to 
the probabilities for flavour transitions Vf, —> i/f m given by 

P{vfi -> v fm ) = Si m - 4^ReJ^ i/m sin 2 -2^IrnJ^ i/m sin2A ij (1) 

i>j i>j 

' -V-' S -V-' 

Pcp Pj2tf 

where the shorthands := UuU*-UmiUmj and Ay := have been used. These generalized vacuum 

transition probabilities depend on all combinations of quadratic mass differences Am 2 , = m 2 — m 2 as well 
as on different products of elements of the leptonic mixing matrix U. We will assume a three neutrino 
framework, i.e. 1 < i,j < 3 and U is a 3 x 3 mixing matrix parameterized in the standard way 

( C12C13 S12C13 si 3 e~ lS \ 

~ S 12 C 23 — Ci2S23Sl3e l<5 C 12 C 23 — Si2S23Si3e l<S S 23 C 13 1 , (2) 

S 12 S 23 — Ci2C23Si3e l<S — C 12 S 23 — Si2C23Sl3e l<5 C 23 C 13 J 

where Cij = cos (Oij) and Sij = sin((%). U contains three leptonic mixing angles and one Dirac-like leptonic 
CP phase <5. Note that the most general mixing matrix for three Majorana neutrinos contains two further 
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Majorana-like CP phases, but it can easily be seen that these extra phases do not enter in the above 
oscillation formulae. Disappearance probabilities, i.e. the transitions Uf t —> Vf i: do not even depend on <5, 
since j( l C l is only a function of the modulus of elements of U. Appearance probabilities, like v e —> are 
therefore the place where leptonic CP violation can be studied. Eq. 0 contains a CP conserving part Pep 
and a CP violating part Pjpp, and both terms depend on the CP phase S. An obvious extraction strategy for 
CP-violation would thus be to look at CP asymmetries m Note, however, that the beams of a long baseline 
experiment traverse the Earth and the presence of matter violates CP by itself. This implies modifications 
of eq. 0 and it makes a measurement of leptonic CP violation more involved. 

The general oscillation formulae in vacuum, eq. 0 , lead to well known, but rather lengthy trigonometric 
expressions for the oscillation probabilities in vacuum. These expressions become even longer and do not 
exist in closed form when arbitrary matter corrections are taken into account. For effectively constant matter 
densities, which is often a good assumption, the problem simplifies somewhat, but the general oscillation 
probabilities are still very lengthy. The Hamiltonian describing three neutrino oscillation in matter can then 
be written in flavour basis as 


(3) 


The first term describes oscillations in vacuum in flavour basis. The quantities A and A' in the second term 
are given by the charged current and neutral current contributions to coherent forward scattering in matter. 
The charged current contribution is given by 


H= — U 1 

( ™i 
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(4) 


where Gf is Fermi’s constant, Y is the number of electrons per nucleon, m n is the nucleon mass and p is the 
matter density. A is positive for neutrinos in matter and anti-neutrinos in anti-matter, while it is negative for 
anti-neutrinos in matter and neutrinos in anti-matter. The flavour universal neutral current contributions 
A! lead to an overall phase which does not enter the transition probabilities. The over-all neutrino mass 
scale rn'i can be written as a term proportional to the unit matrix and can similarly be removed, such that 
only A rri^i and Am§! remain in the first term of eq. ©• After re-diagonalization of the Hamiltonian in 
constant matter density one finds that matter effects lead in a very good approximation to an A-dependent 
parameter mapping in the 1-3 subspace which can be written as m 

sin 2 26 »i 3 


sin 20 i 3 ;m — 


c 2 ± 


A ™ll ,m 

Am 32 ,m 


ATO91 


A m\ x C± , 

Am§! (C± + 1) + A 
2 

Am§! (C± - 1) - A 


(5) 

( 6 ) 

(7) 

( 8 ) 


The index m denotes effective quantities in matter, where 


C± = 


f A 


\ Am| 


31 


— cos 28 ) + sin 2 28 . 


(9) 
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Note that A in C± can change its sign and the mappings for neutrinos and anti-neutrinos are therefore 
different, resulting in different effective mixings and masses. This is an important effect, which will allow 
detailed tests of coherent forward scattering of neutrinos in matter. Note that oscillations in matter depend 
unlike vacuum oscillations on the sign of Am^. This allows to extract the sign(Am'^ 1 ) via matter effects. 

Inserting the parameter mappings eqs. ©-10 into the full oscillation formulae leads still to quite lengthy 
expressions for the oscillation probabilities in matter, where it is not easy to oversee all effects. It is therefore 
instructive to simplify the problem further such that a qualitative analytic understanding of all effects be¬ 
comes possible, while quantitative statements should be evaluated numerically with the full expressions. The 
key for further simplification is to expand the oscillation probabilities in small quantities. These expansion 
parameters are a = Am^/Am^ =O(10“ 2 ) and sin 2 2613 < 0.16. The matter effects can be parameterized 
by the dimensionless quantity A = A/Am^ = 2VE/Am^, where V = \j2GFn & . The oscillation probabil¬ 
ities for all channels can in this way be significantly simplified ESI- Using A = A 31 , the leading terms for 
P{vy. —> Up) and P{y e —> v,j) can, for example, be written as fT711231122| 

P{y^ -> V,/) « 

1 — cos 2 9 13 sin 2 2023 sin 2 A + 2 a cos 2 633 cos 2 0 i2 sin 2 20 23 A cos A, (10) 


P(v e —> v^) « sin 2 2013 sin 2 0 23 


sin 2 (( 1 — A) A) 


± sin S ■ sin 20 13 a sin 20 i2 cos 0 i 3 sin 20 23 sin( A) 

+ cos 6 ■ sin 2013 ce sin 20 i2 cos 0 i 3 sin 20 23 cos(A) 
sin 2 (HA) 


sin(AA) sin((l— A) A) 
A(l-A) 

sin(HA) sin((l— H)A) 

A0--A) 


+ a sin 20 i 2 cos 0 2 3 - 


A 2 


( 11 ) 


where “+” in eq. m stands for neutrinos and for anti-neutrinos. The most important feature of 
eq. dill) is that all interesting effects in the v e —* transition depend crucially on 0 i 3 . The size of sin 2 2013 
determines thus if the total transition rate, matter effects, effects due to the sign of Aand CP violating 
effects are measurable. This is the reason why the size of 0 i 3 is one of the most important questions for 
future oscillation experiments. 

Before we discuss some features of eqs. m and m in more detail, we would like to comment on the 
underlying assumptions and the reliability of these equations. First eqs. El and CD are an expansion in 
terms of the small quantities a and sin 20 i 3 . Higher order terms are suppressed at least by another power 
of one of these small parameters and these corrections are thus typically at the percent level. Note that the 
expansion in a is actually an expansion in the solar and not the atmospheric frequency. The expansion does 
therefore not break down at the first atmospheric oscillation maximum, i.e. at A ~ 1, but at much larger 
baselines before the first (sub-dominant) solar oscillation maximum, i.e. at qA ~ 1. The latter condition 
gives an upper bound for the baseline where eqs. El and CD are good approximations 


L < 8000 km [ 


\GeU/ 


10 ~ 4 eV 2 

Am 2 ! 


( 12 ) 
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while the first oscillation maximum sits at a-L ~ L/30. Eqs. (11 011 and CD are therefore excellent approxima¬ 
tions at and well beyond the first oscillation maximum of long baseline experiments. The matter corrections 
in eqs. 113 and CD are derived for constant average matter density which is a good approximation. 

Note that all quantitative results which will be presented are based on numerical simulations in matter. 
The results are therefore not affected by any approximation. Eqs. 113 and CD will only be used to 
understand the problem analytically, which is extremely helpful in order to oversee the multi-dimensional 
parameter space. 

In addition to long baseline experiments, reactor experiments with identical near and far detectors have 
an excellent potential for precise measurements. The near detector is used to eliminate many common 
systematical errors and the far detector is located typically at a baseline of a few kilometer. For these short 
baselines matter effects can be ignored and one finds to second order in the small quantities sin 29i 3 and a 
for the oscillation probability 

1 — P S e = sin 2 2 0i 3 sin 2 A 31 + a 2 A 31 cos 4 9 i3 sin 2 20 12 . (13) 

At the first atmospheric oscillation maximum, A 31 is approximately tt/2 and sin 2 A 31 is close to one, which 
means that the second term on the right-hand side of this equation can be neglected for sin 2 29i 3 > 10 ” 3 . 
The reactor measurement is dominated in this case at short baselines by the product of sin 2 29 33 and sin 2 A 31 , 
which must be measured as deviation from one. Eq. CD implies that correlations and degeneracies play 
essentially no role in reactor experiments. The behavior in the sin 2 20 13 -A? 7 i 2 1 -plane will also be different 
since eq. C3 is essentially independent of AA reactor experiment will improve the global parameter 
determination in two ways. First, a direct, essentially uncorrelated and clean measurement for 9 33 1251 can 
be obtained which can be used to disentangle the long baseline results. Secondly, the reactor measurement 
can replace the cross-section suppressed anti-neutrino running of the accelerator experiments, leading to 
statistical improvements in the neutrino measurements plfil . 

3 Correlations and degeneracies 

Eqs. 1101 and ( 111 ! exhibit certain parameter correlations and degeneracies, which play an important role 
in the analysis of long baseline experiments, and which would be hard to understand in a purely numerical 
analysis of the high dimensional parameter space. The most important properties are: 

• Eqs. itTTl and (fTTl depend only on the product a • sin 29\ 3 or equivalently Arn 21 • sin 29\ 3 . This are the 
parameters related to solar oscillations which will be taken as external input. Note that the product 
is better determined than the product of the individual measurements of Am 21 and sin 20 i 2 . 

• Next we observe in eq. CD that the second and third term contain both a factor sin(AA), while the 
last term contains a factor sin 2 (A A). Since A A = VL/ 2, we find that these factors depend only on 
L , resulting in a “magic baseline” m when VL mag i C = 27r, where sin(AA) vanishes. At this magic 
baseline only the first term in eq. CD survives and P{v e —> v^) does no longer depend on S , a and 
sin 2012 - This is in principle very important, since it implies that sin 2 29\ 3 can be determined at the 
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magic baseline from the first term of eq. CD whatever the values and errors of 5, a and sin 2 0 12 are. 
For the given matter density of the Earth we find L mag i C = 2tt/V — 8100 km which fits nicely into 
the Earth. This value is quite amazing, since V is given in terms of completely unrelated constants of 
nature like Gf- 

• Next we observe that only the second and third term of eq. CD depend on the CP phase S, and both 
terms contain a factor sin 2013 • a, while the first and fourth term of eq. CD do not depend on the 
CP phase S and contain factors of sin 2 29±3 and a 2 , respectively. The extraction of CP violation is 
thus always suppressed by the product sin 2013 • a and the CP violating terms are obscured by large 
CP independent terms if either sin 2 2013 <C a 2 or sin 2 2013 3> a 2 . The relative contribution of the CP 
phase S to the probability is thus largest for sin 2 20i 3 ~ 40 2 3 ~ a 2 . 

• Another observation is that the last term in eq. CD- which is proportional to a 2 = (A?n| 1 ) 2 /(A?n 2 1 ) 2 , 
dominates in the limit of tiny sin 2 20 i 3 . The error of Am^ limits therefore for small sin 2 29 13 the 
parameter extraction. This last term implies a finite transition probability even for $13 = 0. Observing 
i/ e —> i/fj, or l/fj, —> v e appearance transitions does therefore not necessarily establish a finite value of 
013 = 0 in a three flavour framework. 

• Eqs. CD and CD suggest that transformations exist which leave these equations invariant. Therefore 
degeneracies, i.e. parameter sets having identical oscillation probabilities for a fixed L/E are expected. 
An example of such an invariance is given by a simultaneous replacement of neutrinos by anti-neutrinos 
and Am 2 ! —> —Am^. This is equivalent to changing the sign of the second term of eq. (Hill and 
replacing a —> — a and A —► — A, while A —> A. It is easy to see that eqs. CD and CD are unchanged, 
but this is not a degeneracy, since neutrinos and anti-neutrinos can be distinguished experimentally. 

• The first real degeneracy can be seen in the disappearance probability eq. CD- which is invariant 

under the replacement 023 —> 7 r /2 — 023- Note that the second and third term in eq. ( 1111 ) are not really 
invariant under this transformation, but this change in the sub-leading appearance probability can 
approximately be compensated by tiny parameter shifts. This implies that the degeneracy can in 
principle be lifted with high precision measurements in the disappearance channels. 

• The second degeneracy can be found in the appearance probability eq. CD in the (5— 0 i 3 )-plane !5D . 
In terms of 0 i 3 (which is small) and <5 the four terms of eq. (Hill have the structure 

P(v e -*■ Vp) « 0 2 3 • Fi + 0 i 3 • (± sin 5F 2 + cos 6 F 3 ) + F 4 , (14) 

where the quantities F t . i = 1, ..,4 contain all the other parameters. The requirement P{y e —> = 

const, leads for both neutrinos and anti-neutrinos to parameter manifolds of degenerate or correlated 
solutions. Having both neutrino and anti-neutrino beams, the two channels can be used independently, 
which is equivalent to considering simultaneously eq. CD for F 2 = 0 and F 3 = 0. The requirement 
that these probabilities are now independently constant, i.e. P(y e —> = const, for F 2 = 0 and 

F3 = 0 , leads to more constraint manifolds in the (5 — 0 i 3 )-plane, but some degeneracies still survive. 
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• The third degeneracy m is given by the fact that a change in sign of Am § 1 can essentially be 
compensated by an offset in S. Therefore we note again that the transformation Ajtj^ —> —Am^ 
leads to a —> —a, A —> —A and A —> —A. All terms of the disappearance probability, eq. m , are 
invariant under this transformation. The first and fourth term in the appearance probability eq. HU, 
which do not depend on the CP phase <5, are also invariant. The second and third term of eq. IB 
depend on the CP phase and change by the transformation Am^ —> — Am^. The fact that these 
changes can be compensated by an offset in the CP phase <5 is the third degeneracy. 

• Altogether there exists an eight-fold degeneracy m, as long as only the v e —> 

and v e —> channels and one fixed L/E are considered. However, eqs. (tToll and IB also imply that 

the degeneracies can be broken by using in a suitable way information from different L/E values. This 
can be achieved in total event rates by changing or combining different L or E EUCHIIB but it can 
in principle also be done by using information in the event rate spectrum of a single baseline L, which 
requires detectors with very good energy resolution [B Another strategy to break the degeneracies 
is to include further oscillation channels in the analysis (“silver channels”) |J3J SSI ■ 

The discussion of this section shows the strength of the analytic approximations, which allow to understand 
the complicated parameter interdependence. It also helps to optimally plan experimental setups and to find 
strategies to resolve the degeneracies. 

4 The potential of future neutrino oscillation experiments 

Triggered by the spectacular results in neutrino physics during the last ten years, several new experimental 
projects are under way in this field. It is therefore interesting to investigate where we should stand in the 
determination of neutrino oscillation parameters in ten years from now. It is also interesting to look further 
and to estimate the ultimate precision which could be obtained. 

The precision of quantities like sin 2 2Qi% which is found form the simulation of experiments will be 
presented in a way shown in fig. ^ The bands show how the initial value, which is given by statistics alone 
(left edge of blue/dark grey band) deteriorates by systematic errors, by parameter correlations (e.g. with 
the unknown or partly known CP phase) and parameter degeneracies (due to trigonometric ambiguities). 
It is important to note that a given experiment (or combination of experiments) typically measures some 
parameter combination with a precision which is considerably better than the final limit. This precision 
of the experiment is shown in fig. |2 as the right edge of the blue/dark grey band. This precision might 
be called (sin 2 2 ^ 13 ) e ^, since it expresses the precision if all other unknown parameters are fixed and no 
errors are included. However, if one properly extracts a limit of sin 2 2013 with all unknowns properly taken 
into account, then one ends up at the right edge of the yellow/light grey band. Distinguishing in this way 
between the precision and the sensitivity is quite useful, since it also shows the room for improvement by 
combinations with other similarly precise experiments with other parameter dependence. 
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statistical limit 
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Figure 1: The precision for sin 2 2013 is shown in colored bands, where the left edge of the blue/dark grey band 
shows the initial value which is obtained if only statistics is considered. The right edge of the blue/dark 
grey band is the result after the systematic errors are included. This is the principal precision of the 
experiment. However, the sensitivity for sin 2 2013 deteriorates further due to parameter correlations and 
parameter degeneracies. The final value is the right edge of the yellow/light grey band. The range covered 
by the green and yellow bands can lead to remarkable synergies when this experiment is combined with 
another experiment with similar precision, but different parameter dependence. 


4.1 Next generation experiments 

Future oscillation experiments can be grouped according to their time scale of operation. The K2K experi¬ 
ment is already running and it tests the leading atmospheric oscillation already now. Next come the MINOS 
and CNGS projects which are already operating or under construction, respectively. Therefore we include 
in our study the conventional beam experiments MINOS UJJ, and the CNGS experiments ICARUS 23 
and OPERA [23. We include also the subsequent superbeam experiments J-PARC to SuperKamiokande 
(T2K) III and NuMI off-axis (NO^A) Rlj . as well as new reactor neutrino experiments (42) with a near 
and far detector. The main characteristics of these experiments are summarized in tab.^ For the reactor 
experiments we use the Double-CHOOZ proposal (D-CHOOZ) m as initial stage setup with roughly 6 x 10 4 
events, and an optimized setup called Reactor-II, with a slightly longer baseline and 6 x 10 5 events. Such 
a configuration could be realized at several other sites under discussion 221- The results presented in the 
following are based on Ref. EH, where more details on the analysis can be found. The simulations of the 
experiments as well as the statistical analysis is performed with the GLoBES software package E3 


Label L [km] 

{Eu) 

^run 

channel 

Conventional beam experiments: 


MINOS 735 

3 GeV 

5yr 

v i*~ 


ICARUS 732 

17 GeV 

5yr 



OPERA 732 

17 GeV 

5yr 

v i*~ 


Off-axis superbeams: 




T2K 295 

0.76 GeV 

5yr 

v i*~ 

^e,fi 

NO^A 812 

2.22 GeV 

5yr 

V H~ 

V e,fi 

Reactor experiments: 




D-CHOOZ 1.05 

~ 4MeV 

3yr 

v e - 

->v e 

Reactor-II 1.70 

~ 4MeV 

5yr 

v e - 



Table 1: Characteristics of the considered experiments. 









A first interesting question concerns improvements of Am § 1 and sin 2 0 2 3 - I n tab. Owe show the precision 
which can be obtained in the future in comparison to the current precision, as obtained from a global fit to 
SuperKamiokande (SK) atmospheric and K2K long-baseline data dSl. The last row is the precision which 
can be obtained by combining all experiments. We observe from these numbers, that the accuracy on A 
can be improved by one order of magnitude, whereas the accuracy on sin 2 0 2 3 will be improved only by a 
factor two. 



|Am 2 3 | 

sin 2 023 

current 

88 % 

79% 

MINOS+CNGS 

26% 

78% 

T2K 

12 % 

46% 

NOz/A 

25% 

86 % 

Combination 

9% 

42% 


Table 2: Precision for |Am 31 | and sin 2 (9 23 at 3 a for the values Am^ = 2 x 10 3 eV 2 , sin 2 0 23 = 0.5. 

Tab. |3 depends on the value of Awhich is shown in fig. O The sensitivity suffers for all experiments 
for low values of Am^. T2K will provide a precise determination of Am^j at the level of a few percent for 
Am|i > 2 x 10 - 3 eV 2 . Although NOi'A can put a comparable lower bound on Am 2 ^ the upper bound 
is significantly weaker, and similar to the bound from MINOS. The reason for this is a strong correlation 
between Am 31 and 0 23) which disappears only for A> 3 x ICO 3 eV 2 . From the right panel of fig. 0 
one can see that for Am 31 ~ 2 x 10 -3 eV 2 only T2K is able to improve the current bound on sin 2 0 2 3 . One 
reason for the rather poor performance on sin 2 62 3 is the fact that these experiments are sensitive mainly to 
sin 2 26 l 2 3. This implies that for 623 ~ 7 r /4 it is very hard to achieve a good accuracy on sin 2 0 2 3 , although 
sin 2 2023 can be measured with relatively high precision ED- 


2 2 

Am 31 -precision sirf0_, 3 -precision 



Figure 2: The precision of Am^ (left panel) and sin 2 623 (right panel) as a function of the “true value” of 
Amf n for 6 % 3 ue = tt/ 4 (from Ref. gZj). 
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Another interesting question is if the next generation long baseline experiments which will operate in 
the next years will be able to test the three flavourness of the oscillations. The sensitivity to a finite value 
of the key parameter $13 is shown in fig. [I] for MINOS, OPERA and ICARUS. It can be seen that these 
experiments have only a modest potential for improvements of the existing $13 limit. 



5 10 15 20 25 30 

Luminosity [10 20 pot] 


Sensitivity to si n 2 2C| 3 at 90% CL 



210 -2 io _1 


sin 2 2013 


Figure 3: Left plot: The sensitivity of the MINOS experiment to 6*13 as a function of the protons on target 
(pot) assuming a 5% flux uncertainty. The dashed lines represent what 1,2 and 5 years of operation might 
achieve (from left to right). Right plot: Comparison of 5 years of operation for the MINOS and CNGS 
experiments. The grey area for large sin 2 indicates in all cases the current limit from the CHOOZ 
experiment. The color code of the error bars is explained in fig. ^ Further details can be found in m 


The combined sensitivity to sin 2 20 13 of the next-to-next generation experiments is compared in the left 
panel of fig. Q] with new reactor experiments, T2K (JPARC-SK) and NO^A (NuMI). It can be seen that 
the sin 2 20 i 3 -limits from beam experiments are strongly affected by parameter correlations and degenera¬ 
cies, whereas reactor experiments provide a “clean” measurement of sin 2 20 i 3 , dominated by statistics and 
systematics EH! 



MINOS& 

ICARUS& 

OPERA 

D-CHOOZ 
Reactor II 


JPARC-SK 


NuMI 


CHOOZ & Solar 
excluded (90% CL) 


1(T 2 KT 1 

sin 2 2013 sensitivity limit 


sin 2 20 i 3 sensitivity limit 



sin 2 20 b sensitivity limit 


Figure 4: Left plot: Sensitivity to sin 2 2013 at 90% CL for Am^ = 2 x 10 3 eV 2 , Am^i = 7 x 10 5 eV 2 . 
Right plot: Sensitivity to sin 2 2013 at 90% CL as a function of the “true value” of Am^. 
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The dependence of the sin 2 2 (?i 3 -limit on the value of Ais shown in the right panel of fig. 01 where the 
sensitivity of all experiments gets again rather poor for low values of A. For Ato^ ~ 2 x 10 ~ 3 eV 2 we find 
roughly an improvement by a factor 2 from conventional beam experiments (MINOS+ICARUS+OPERA 
combined), a factor 4 from D-CHOOZ, and a factor 6 from the superbeams T2K and NOizA with respect to 
the current bound from global data m, Note that an optimized reactor experiment such as Reactor-II has 
the potential for even better sin 2 20 i 3 -sensitivities than the superbeams (c.f. left panel of fig. 0 ). 

4.2 Synergies 

The previous discussion shows that competing plans with similar potential might be realized at the same time 
scale. This allows to combine the statistics of similar experiments leading to improved global fits. However, 
it is also possible to utilize synergies between experiments which are more than the simple addition of 
statistics. The point is that individual experiments measure a certain parameter combinations which contain 
different degeneracies and correlations. Experiments with similar sensitivities and different correlations and 
degeneracies allow to separate the parameters partly or fully. An example of such a discussion is given by 
combining the T2K and NOjzA experiments for a fixed time of operation in the best possible way. T2K 
is essentially insensitive to matter effects, while matter effects play already some role for the longer NOizA 
baselines. Both experiments could run partly with neutrino and partly with anti-neutrino beams. The 
cross-sections for anti-neutrinos are, however, smaller, leading to fewer events for the same running period. 
An anti-neutrino running is moreover in many aspects like a different experiment, but it is clear that anti¬ 
neutrino information is crucial in order to resolve the parameters. A comparable reactor neutrino experiment 
would be very useful here. It could provide the required information such that both T2K and NOizA could 
initially run fully in neutrino mode. 

Such a synergetic combination would be especially interesting if sin 2 20\3 would be close to the current 
bound. In order to demonstrate these synergy effects we assume that sin 2 29\3 = 0.1 and investigate what 
we could be learned about the CP-phase S and the neutrino mass ordering. First we note that T2K, NOizA 
and the reactor experiment will all be able to establish the non-zero value of sin 2 29±3. However, depending 
on the unknown value of <5 different values of sin 2 29±3 will be allowed. This can be seen as allowed regions 
in the 0 i 3 -< 5 -plane shown in Figs. 8 and 9 of Ref. m- None of the experiments on their own can give 
any information on the CP-phase S and on the mass hierarchy. The determination of sin 2 26*13 from beam 
experiments is strongly affected by the correlations with 5, and especially for NOtzA also correlations with 
other parameters are important. Moreover, the inability to rule out the wrong mass hierarchy leads to a 
further ambiguity in the determination of sin 2 26 * 13 . In contrast, since the f/ e -survival probability does not 
depend on 6, Reactor-II provides a clean determination of sin 2 26*13 at the level of 20% at 90% CL. If all 
experiments are combined the complementarity of reactor and beam experiments allows to exclude up to 
40% of all possible values of the CP-phase for a given hierarchy. The wrong hierarchy can be ruled out at a 
modest confidence level with Ay ; 2 — 3 due to matter effects in NOizA. However, at high confidence levels all 
values of S are allowed, and moreover, even for a given hierarchy CP-conserving and CP-violating values of 
S cannot be distinguished at 90% CL. These results depend also to some extent on the value of 5. 


11 


So far we have considered only neutrino running for the superbeams, since it is unlikely that significant 
data can be collected with anti-neutrinos within ten years from now. Nevertheless, it might be interesting 
to investigate the potential of a neutrino-antineutrino comparison. In fig. [3 we show the results from 
T2K+NOi'A with 3 years of neutrinos plus 3 years of anti-neutrinos each (left), in comparison with the 
case where the antineutrino running is replaced by Reactor-II (right). We find that antineutrino data at 
that level does neither solve the problems related to the CP-phase nor to the hierarchy. Still CP-violating 
and CP-conserving values cannot be distinguished at 90% CL. Moreover, the determination of sin 2 2f?i3 is 
less precise than from the reactor measurement. To benefit from antineutrino measurements a significantly 
longer measurement period would be necessary, to obtain large enough data samples. 


T2K + NOVA T2K + NOvA + Reactor-II 

3 yrs neutrinos + 3 yrs anti-neutrinos T2K + NOvA: 3 yrs neutrinos 



Figure 5: Antineutrino running vs Reactor-II. We show the 90% CL (solid curves) and 3er (dashed curves) 
allowed regions in the sin 2 2di3-<5-plane for the assumed values sin 2 26^ = 0.1 and S = 7t/ 2. The blue/dark 
curves refer to the allowed regions for the normal mass hierarchy, whereas the red/light curves refer to the 
sgn(Am§ 1 )-degenerate solution (inverted hierarchy), where the projections of the minima onto the sin 2 29i3-5- 
plane are shown as diamonds (normal hierarchy) and dots (inverted hierarchy). For the latter, the Ay 2 -value 
with respect to the best-fit point is also given. 


4.3 Long term perspectives 

Beyond the discussed accelerator and reactor based oscillation experiments exist more ambitious projects like 
the JHF-HyperKamiokande project, beta beams and neutrino factories. Such experiments clearly require 
further R&D before they can be built. However, assuming current knowledge, we believe that such setups 
are possible in the long run. The potential of the JHF-HyperKamiokande experiment and a neutrino factory 
are compared in fig. El to T2K (JHF-SK) and NO^A (NuMI). It can be seen that the existing limits can be 
improved by a few orders of magnitude compared to now. 

5 Theoretical implications and conclusions 

One of the most interesting unsolved topics is the origin of flavour and fermion masses. There exist apparent 
regularities in the fermionic field content which make it very tempting to introduce right-handed neutrino 
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Figure 6 : Left plot: The @13 sensitivity of different future accelerator based neutrino oscillation experiments 
EDI • Right plot: The $13 values for which sensitivity to matter effects, i.e. sign{Arn${) exists. The shown 
bands are again the reduction of sensitivity from a purely statistical limit (left end of the dark grey/blue 
range) by systematics (right end of dark grey/blue), correlations (medium grey/green) and degeneracies 
(light grey/yellow). The right end of the light grey(yellow) band represents the final 90%CL limit. The grey 
area for large sin 2 26^ indicates the current limit from the CHOOZ experiment. 


fields leading to both Dirac and Majorana mass terms for neutrinos. Diagonalization of the resulting mass 
matrices yields Majorana mass eigenstates and due to the see-saw mechanism very small neutrino masses. 
This can also be nicely realized in embeddings of the SM into GUTs with larger symmetries, such as SO(IO). 
Before the discovery of large leptonic mixing, many theorists expected the leptonic mixings to be similar to 
quark mixing, characterized by small mixing angles. Experiment led theory in showing the striking results 
that sin 2 2023 — 1 and tan 2 6*12 ~ 0.39, while $13 is small. By finding two large mixing angles, neutrino 
physics has already provided surprising and very valuable information which severely constrains models of 
neutrino masses. Future precision neutrino oscillation experiments will provide further precision tests of the 
flavour sector. The level of precision will confirm or rule out ideas about the origin of flavour and connected 
topics. 

An important subject is the small value of sin 2 20 i 3 . Since there are two large mixing angles, there is 
no particular reason to expect the third angle, $ 13 , to be extremely small or even zero. A small value of 
sin 2 2013 — 0.1 could be a numerical coincidence in a framework which predicts generically large or sizable 
mixings. However, if the limit on sin 2 2013 would become smaller by an order of magnitude then some 
protective mechanism like a “symmetry” would be required. This can be seen in neutrino mass models 
which are able to predict the large values for 0 3 2 and 023- Such models have a certain tendency to predict a 
sizable value of 0 i 3 as can be seen, for example, in g3EniiniEllE3E2 The conclusion is that a value of 
013 close to the CHOOZ bound would be quite natural, while much smaller values are less likely or hard to 
understand. 

Future precision measurements can also test if relations like 023 = 7 t /4 ^5J or 0 i 2 + 0c = 7 r /4 which 
are currently fulfilled within experimental errors still hold at much better precision. If so, then this would 
provide strong constrains on the origin of the flavour structure. Precision is also valuable, even if such special 
relations are not found. The point is that it is generally not easy to predict a set of very precisely known 
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masses and mixings in a certain model or class of models. 

Neutrino masses and mixing parameters are also subject to quantum corrections between low scales, 
where measurements are performed, and high scales where some theory typically predicts the masses and 
mixings. This has interesting implications, since it implies that certain deviations from special relations are 
expected due to quantum corrections (renormalization group or RGE effects). Suppose, for example, that 
some theory were able to predict 6*13 = 0. Then RGE effects would still predict a tiny, but finite value at low 
energy. Strictly speaking, $13 = 0 cannot be excluded completely by this argument, as the high-energy value 
could be just as large as the change due to running and of opposite sign. However, a complete cancellation 
of this kind would be a miraculous fine-tuning, since the physics generating the value at high energy is not 
directly related to quantum corrections at lower energies. The strength of the running of $13 depends on the 
neutrino mass spectrum and whether or not supersymmetry is realized. For the Minimal Supersymmetric 
Standard Model one finds a shift A sin 2 29\^ > 0.01 for a considerable parameter range, i.e. one would 
expect to measure a finite value of $13 1561 . Conversely, limits on model parameters would be obtained if 
an experiment were to set an upper bound on sin 2 2 f?i 3 in the range of 0.01. In any case, it should be clear 
that a precision of the order of quantum corrections to neutrino masses and mixings is very interesting in a 
number of ways. 

Precision measurements would also allow very interesting tests of many other topics, like MSW matter 
effects, three neutrino unitarity, neutrino decay, de-coherence and NSI effects. There exists also an inter¬ 
esting interplay with theories beyond the Standard Models, flavour models as well as astroparticle physics 
(leptogenesis, supernovae, nucleosynthesis, structure formation). In summary, future precision oscillation 
experiments will provide the best window into the so far un-understood flavour sector. It may give us a 
glimpse on the origin of flavour, but it may also lead to unexpected results and insights as it happened before 
in neutrino physics. 
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